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EXECUTIVE  SUMMARY 

Rotating  machinery  comprises  the  most  important  class  of  systems  associated  with  power  plants,  transportation, 
and  manufacturing.  With  the  dramatic  increase  in  technology  and  economic  pressure,  the  performance  demands 
related  to  rotating  machinery  have  made  it  essential  that  advanced  diagnostic  systems  capable  of  early  detection 
and  identification  of  faults  and  changing  conditions  associated  with  reduced  performance  be  developed.  These 
diagnostic  systems  must  not  only  take  advantage  of  the  most  advanced  theories,  but  they  must  be  user-friendly. 
The  value  of  advanced  techniques  is  all  too  often  compromised  by  (i)  failure  to  take  real  world  conditions  which 
might  violate  theoretical  assumptions  into  account  and  (ii)  the  inability  of  applications  engineers  to  apply  them 
to  real  world  problems.  Having  said  this,  the  goal  of  the  research  program,  the  results  of  which  are  summarized 
in  this  report,  was  threefold:  first,  to  dev  elop  a  rigorous  mathematical  approach  that  would  take  greater 
advantge  of  the  quasi-periodic  structure  of  phenomena  such  as  rotating  machinery;  second,  to  integrate  the 
development  of  this  approach  with  particular  real  world  applications  which  might  serve  to  suggest  the  extent  to 
which  assumptions  related  to  the  theory  being  developed  might  hold  in  practice,  and  hence  guide  the  research 
direction;  and  third,  to  develop  research  collaborations  with  researchers  in  other  disciplines  who  might  be  able 
to  provide  new  and  valuable  perspectives,  and  who  might  also  directly  value  from  the  development.  Not  only 
does  this  collaboration  tend  to  preclude  the  “re-invention  of  the  wheel"  at  a  time  when  it  is  extremely  difficult 
for  researchers  in  one  area  to  keep  apprised  of  related  mathematical  problems  and  results  in  other  disciplines,  it 
also  supports  cross  fertilization  which  can  be  mutually  beneficial 

The  approach  that  was  taken  in  this  effort  was  based  on  the  fact  that  very  often  periodic  stationary  processes 
include  two  distinct  components,  namely,  point  spectrum  associated  with  deterministic  processes,  and 
continuous  spectrum  associated  with  regular  random  processes.  Interestingly,  traditional  methods  of  spectral 
analysis,  which  have  been  used  in  relation  to  rotating  machinery  for  decades,  are  in  fact  ill-conditioned  in 
relation  to  random  processes  having  such  a  mixed  spectrum.  Another  apsect  of  our  approach  which 
distinguishes  it  from  not  only  traditional  methods,  fouA  foora  all  methods  proposed  to  date,  is  the  use  of  an  entire 
family  of  spectral  estimators,  as  opposed  to  a  single  best  estimator,  to  characterize  the  spectral  information.  The 
mathematical  foundation  for  this  approach  was  developed  prior  to  this  research  effort.  The  major  results  of  this 
research  effort  include  the  following: 

(i)  characterization  of  the  influence  of  period  uncertainty  on  estimation  of  a  periodic  time/frequency  spectrum 
associated  with  a  nominally  wide  sense  cyclostationary  (wsc)  process, 

(ii)  large  sample  distribution  descriptions  for  the  AR(p)  and  M V(p)  spectral  estimators  for  processes  with  mixed 
spectrum. 


(iii)  a  time-to-angle  transformation  method  to  better  accommodate  period  variability,  combined  with  an 
improved  method  for  tracking  real  sinusoids  with  slowly  varying  frequency. 

(iv) greater  insight  into  issues  related  to  application  of  advanced  spectral  analysis  methods  for  characterizing 
random  processes  associated  with  engines,  compressors,  and  helicopter  drivetrains, 

(v) a  Matlab-based  virtual  signal  analyzer  which  incorporates  a  number  of  our  results  in  a  user-friendly  fashion, 

(vi)  development  of  research  collaborations  and  workshops  w  hich  serve  to  bring  signal  processing  problems 
associated  with  rotating  machinery  to  a  broader  base  of  researchers  in  industry,  defense  and  academic 
institutions. 
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€.  SUMMARY  OF  RESEARCH  RESULTS 
Cl.  BACKGROUND 

A  discrete  time  random  process,  Xt,  is  said  to  be  a  mixed  process  if  Xt  —  St  +  Nt,  where  St  is  a  sum  of 
sinusoids  and  where  Nt  has  t  +  r)  which  is  square  summable  in  r  for  any  fixed  t.  In  the  case  where  Nt  is 
wsj,  define  the  unbiased  lagged-product  autocorrelation  estimator 


Rx{r)  =  (1/NYE  XtXt+r.  (1) 

t=0 

The  following  large  sample  result  obtained  in  1994  by  Li  et  al  [LKY]  for  a  mixed  wss  process,  Xt  is  the  basis 
of  our  key  theoretical  results. 

Theorem  1.  [LKY].  Let  xt  be  a  wss  process  with  point  spectrum  at  and  with  a  Gaussian  regular  process. 
Let  rp  —  [r( 0), r(p)]  with  r(k)  in  (2),  and  let  7P  be  the  unbiased  lagged-product  estimator  of  rp,  Then 
yJN{  rp  -  rp )  ~  N( 0,  E)  where  E  =  {crjtk} 

q  oo  oo 

where  crjtk  =  2A*lcos(jLJm)cas(ku>m)]J2  r(r)cos(u>mT)  +  ]T  iTiT)T(T  +  j-k)  +  r(r  +  j) r(r  -  k)}. 

771=1  T—  CO  r=— 00 

The  above  theorem  is,  in  out  belief,  a  major  breakthrough,  in  that  it  is  the  first  such  result  related  to  the  use  of 
the  lagged-product  autocorrelation  estimator  for  mixed  processes.  It  is  well  known  that  sinusoids  appear  as 
Dirac  delta  functions  in  a  spectral  density  representation.  Practically  speaking,  they  contribute  to  a  spectral 
density  approximation,  based  on  a  finite  number  of  correlation  lags,  peaks  whose  widths  and  amplitudes  are  a 
function  of  the  number  of  lags  used.  In  fact,  as  the  number  of  lags  approaches  infinity,  the  peaks  become 
unbounded.  Consequently,  if  one  is  unaware  of  the  presence  of  sinusoids,  then  one  may  be  easily  misled  in 
conducting  spectral  analysis  and  system  identification  (see,  e.g.  [LSLj).  It  follows  that  an  essential  ingredient  of 
analysis  of  any  spectral  analysis  methodology  for  processes  which  have  mixed  spectrum  is  a  method  to 
decompose  the  process  being  studied  into  its  continuous  and  point  spectrum  components.  A  limited  number  of 
methods  have  been  proposed  for  this  purpose,  the  majority  of  which  assume  the  regular  component  to  be  white. 
The  method  of  Thompson  [THO],  which  is  perhaps  the  most  well  known,  does  not  require  the  regular  portion  to 
be  white.  The  convergence-based  method  of  [FFS],  w'hich  we  have  found  to  be  valuable  in  a  wide  variety  of  real 
world  settings,  is  the  method  we  have  focused  on  in  this  work  for  detecting  sinusoids.  It  has  a  number  of 
advantages  over  Thompson's  method.  One  is  that  it  does  not  rely  on  a  single  spectral  estimate,  but  rather  on  an 
entire  family  of  estimates.  This  allows  for  a  much  more  visual  interpretation.  Second,  and  more  important,  it  is 
equally  valid  for  a  large  class  of  nonstationary  processes  [FRS].  The  family  of  spectra  studied  in  [FFS]  have 
been  labeled  variously,  the  Capon,  maximum  likelihood,  and  minimum  variance  (MV)  spectra.  Here,  we  will 
use  the  MV  terminology.  The  MV(rc)  spectrum  for  a  scalar  process  Xt  is  defined  as 

MVn(u>)  4  (e-i*-^)-1  =  (2) 

0 

where,  e=e(w)=[l,  e^, ...,  einu,]tr,  R  —  {R(j  -  k);fk  =  1, ...,  n},  and  *  denotes  the  conugate  transpose.  It  is 
attractive  from  a  spectral  analysis  perspective,  since  it  implies  that  AR  and  MV  spectra  may  be  computed 
jointly,  with  little  additional  computational  effort  When  first  proposed  by  Capon  [CAP],  and  subsequently 
shown  by  Lacoss  [LAC]  to  have  less  resolving  power  than  the  AR  spectrum,  the  MV  spectrum  was  generally 
viewed  as  a  possible  alternative  to  the  AR  spectrum.  However,  the  result  of  [FFS],  namely  that  (2)  converges  to 
the  point  spectrum  as  n— kx>,  suggests  that  the  AR  and  MV  spectra  actually  complement  each  other  in  mixed 
spectrum  analysis. 
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C2.  THEORETICAL  RESULTS 

The  majority  of  the  results  summarized  in  this  section  were  first  documented  in  [LAU]  and  [LIU].  They  have 
also  been  submitted  for  publication  [LAS],  [LSW],  [LIS].  We  first  present  some  large  sample  results  pertaining 
to  AR  and  MV  spectral  estimators  for  mixed  processes. 

Theorem  2.  [LAU]  Let  Sp(u )  be  the  AR(p)  spectral  estimator  for  xt7  and  let  SP(ui)  be  the  corresponding  AR(p) 
approximant  obtained  from  the  theoretical  Yule-Walker  equations .  Then  ^/N  (Sp(u)  —  Sp(lj))  ~  N(0,^tp). 

The  covariance  expression  is  quite  complicated,  and  so  is  omitted  here  for  convenience.  The  above  result  reverts 
back  to  the  well  known  result  for  a  regular  process  (e.g.  [BRD]),  at  frequencies  sufficiently  removed  from  the 
point  spectrum  frequencies.  At  the  point  spectrum  frequencies  the  mean  becomes  unbounded  as  the  AR  order 
approaches  infinity.  Furthermore,  the  variance  at  these  frequencies  is  not  as  well  behaved  as  at  other 
frequencies.  While  we  have  not  been  able  to  arrive  at  results  as  both  the  model  order,n,  and  data  length,  N 
approach  infinity,  our  simulations  suggest  that  at  the  point  spectrum  frequencies  the  condition  that  n*  /N-> 0,  as 
opposed  to  n2/N->07  is  required  for  the  variance  to  go  to  zero.  A  somewhat  similar  distinction  of  behavior  at 
point  spectrum  frequencies  was  discovered  for  the  MV  spectral  estimator. 

Theorem  3.  [LIU]  Let  MVn (u)  be  the  lagged-product  estimator  ofMVn{cj)  defined  in  (9).  Then 
y/W{MVn{uj)  -MVn(u>)  ~  N(0,Qn(<jj))  where  Qn(u)  =  MV* (u) f3n (u/)trE  & (a;).  The  pxl  vector  (3n(uj)  has  ith 
component  AM  =  e*(u;)R“1Anti_i/?”1enM  where  An,ihas  ±ith  off-diagonals  equal  to  one ,  and  is  zero 
elsewhere. 

The  complex  form  of  the  variance  of  the  MV(n)  estimator  has  thus  far  precluded  insight  into  its  behavior.  For 
this  reason,  we  have  explored  this  behavior  via  a  variety  of  simulations  in  [LIU].  Foremost  among  our 
observations  to  date  is  that  at  non-point  spectrum  frequencies  this  variance  goes  to  zero  as  n-+oo  in  the  same 
manner  as  the  mean  does.  While  it  is  known  that  the  mean  has  a  -3  dB  per  doubling  of  n  rate,  the  variance  rate 
may  not  be  the  same.  At  point  spectrum  frequencies  this  variance  converges  to  a  nonzero  value  which  is 
dependent  on  the  local  signal-to-noise  ratio  (SNR).  While  we  have  proposed  detection  schemes  based  on  this 
family  for  the  random  process  setting  [SHL,  [LYS],  they  are  non-statistical  in  nature.  The  challenge  here  is  the 
same  as  for  any  detection  scheme.  The  more  correlations  one  can  use,  the  better  the  detection  performance.  But 
the  statistical  variability  of  higher  order  lags  also  increases  with  n.  The  above  theorems  provide  the  tools  for 
balancing  these  factors. 

One  of  the  original  goals  of  the  research  was  to  characterize  the  statistics  of  the  AR(p)  estimator  of  the  periodic 
time/frequency  spectrum  associated  with  a  wcs  process.  However,  we  discovered  that  any  amount  of  period 
uncertainty  in  such  a  process  will  have  the  effect  of  destroying  the  time/frequency  information.  This  behavior  is 
summarized  in  the  following  theorem. 

Theorem  4  [SHW].  Let  x(t)  be  a  wcs  process  with  period  Ta,  and  let  S(uj1  t)  be  its  associated  spectrum.  Let 
S(u,  t)  be  an  estimator  of  S(w,t)  based  on  the  autocorrelation  estimator 

^  m 

R(t ,  t  +  r)  =  (1  /m)J2x(kT  +  t)x(kT  +  t  +  r).  Then  if  T  ^  T0  or  if  TQ  has  any  randomness  over  time ,  the 

k=l 

estimator  S(u1  t)  will  become  independent  of  t  as  m  approaches  infinity. 

It  is  the  behavior  described  in  this  theorem  which  guided  a  large  part  of  our  research,  in  an  effort  to  arrive  at 
advanced  spectral  methods  for  processes  such  as  those  associated  with  rotating  machinery  operating  at  a 
nominally  constant  speed,  but  having  a  measurable  amount  of  period  variability.  The  following  theorem  is 
related  to  our  attempt  to  accommodate  it. 
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Theorem  5  [LAS],  Let  x{t)  =  sin(ut  +  6)  +  i /(£),  where  z/(t)  /s  Wz/te.  D  be  the  least  squares  estimator  of  u>, 
based  on  {x( t)}^.  Then  for  large  n,  Q  is  normally  distributed  with  mean  and  variance  described  in  [LAS]. 

The  expressions  for  the  mean  and  variance  are  a  bit  involved,  and  so  are  not  given  here.  The  reader  is  referred 
to  [LAS].  This  theorem  is  significant  for  two  reasons.  First,  it  provides  a  quantitative  rationale  for  positioning  a 
sinusoid  (or  any  narrowband  process)  in  the  center  of  the  frequency  analysis  range.  Such  positioning  has  long 
been  noted  to  yield  better  AR  descriptions  of  the  spectral  information,  but  only  in  a  qualitative  sense. 
Specifically,  not  only  are  both  bias  and  variance  minimized  at  this  position,  but  the  bias  is  exactly  zero.  Second, 
the  theorem  provided  a  guideline  for  development  of  an  improved  time-to-angle  transformation  described  in 
[LAS].  Such  a  transformation  can  not  only  obviate  the  need  for  expensive  hardware  to  accomplish  the  same 
thing,  but  can  be  applied  to  data  that  has  been  praiously  collected  at  uniform  time  sampling  interval.  Finally, 
we  present  a  conjecture  based  on  our  latest  investigation  of  the  properties  of  the  AR(n)  spectral  estimator 
applied  to  a  mixed  process. 

Conjecture  [LS2].  The  AR(p)  spectral  estimator  for  a  mixed  process  converges  to  the  continuous  spectrum 
except  at  frequencies  associated  with  the  point  spectrum ,  where  it  diverges  at  an  asymptotic  rate  of p2. 

This  conjecture  was  arrived  at  after  extensive  simulations.  It  is  important  since  it  is  commonly  assumed  that  one 
can  arrive  at  a  point  spectral  estimate  by  dividing  a  continuous  spectral  estimate  by  p.  This  is  true  of  both  the 
MV(p)  and  p-windowed  autocorrelation  methods.  The  above  theorem  states  that  in  the  case  of  the  AR (p) 
spectrum  one  must  divide  by  p2.  The  practical  implication  of  this  in  relation  to  detection  is  immediate.  By 
dividing  by  p2  one  not  only  converges  to  a  constant  proportional  to  the  SNR,  but  one  reduces  the  continuous 
spectrum  asssociated  with  noise  by  6  dB  per  doubling  of  p,  as  opposed  to  the  traditional  3  dB. 

C3.  APPLICATIONS 

The  following  application  studies  were  conducted  as  part  of  the  research  effort.  They  were  valuable  in 
identification  of  real  world  issues  that  must  be  addressed  to  realize  the  potential  of  advanced  spectral 
techniques. 

Diesel  Engine  Vibration -  Analysis  of  vibration  associated  with  a  diesel  engine  operating  at  nominally  constant 
speed  is  the  subject  of  [SHW].  That  study  led  to  Theorem  4  above,  and  to  other  investigations.  For  example,  it 
pointed  to  some  real  world  complications  concerning  the  information  in  the  MV  and  AR  spectral  estimators.  It 
motivated  an  investigation  of  tone  tracking  algorithms,  which  led  to  Theorem  5.  Finally,  it  demonstrated  the 
difficulty  of  obtaining  useful  time/frequency  spectral  information  in  any  setting  involving  sinusoids  and  period 
uncertainty. 

Compressor  Pressure -  Analysis  of  the  downstream  pressure  associated  with  a  high  speed  compressor  [VFS], 
[FSS]  demonstrated  the  value  of  synchronous  data  collection  in  constructing  time/frequency  spectral 
information.  But  it  also  identified  limitations  of  such  hardware-based  systems,  in  that  they  can  lead  to 
increasing  phase  errors  throughout  the  period.  This  motivated,  in  part,  our  development  of  a  software-based 
time-to-angle  transformation  method  which  can  overcome  this  limitation. 

Helicopter  Gear  Train  Vibration-  Analysis  of  this  data  set.  which  is  provided  by  the  U.S.  Navy  with  the 
assistance  of  Westland  and  Penn  State  University,  provided  us  with  the  opportunity  to  evaluate  our  time-to- 
angle  transformation  method,  and  to  investigate  the  use  of  the  above  results  in  detecting  faults  in  a  well 
controlled  experimental  setting. 

Upper  Atmosphere-  In  an  attempt  to  investigate  the  utility  of  the  convergence-based  MV  method  of  analysis  to  a 
non-engineering  problem,  we  addressed  the  problem  of  modeling  the  upper  atmospheric  wind  profiles  in  a 
region  of  the  South  Pacific.  In  [WIS]  it  is  shown  that  this  approach  led  to  not  only  quantifying  the  amount  of 
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F.  COLLABORATIONS 


FI.  Research  at  Host  Institutions 

Host  Location 


Dates 


1.  Aus.  Natl.  Univ. 

2.  DSTO  (defense  dept) 

3.  DSTO 

4.  DSTO/ANU 

5.  Korea  Adv.  Inst, 
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Adelaide,  Aus. 
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Science  Town 
Taejon,  Korea 
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8/3/92-8/20/92 

7/21/93-8/13/93 

7/95-12/95 

1/19-1/26/96 


F2.  Conference  and  Workshop  Involvement 

•  Invited  to  be  part  of  a  committee  to  organize  a  U.S./Australia  signal  processing  workshop,  held  25- 
27  June  1997.  The  workshop  theme  was  development  and  application  of  adv  anced  signal  processing 
for  industry  and  Defense  applications. 

•  Invited  to  participate  ip  NSF/ONR  workshop  “Signal  processing  for  manufacturing  and  machine 
maintenance",  Alexandria,  VA,  13-15  March  1996. 

•  Invited  to  organize  and  chair  a  special  session  on  rotating  machinery  at  the  next  IEEE  Workshop  on 
Statistical  Signal  and  Array  Processing,  to  be  held  in  Portland,  OR  in  June  1 998. 

•  Principal  organizer  of  an  ONR/AFOSR-sponsored  workshop  on  helicopter  diagnostics  to  be  held 
March  26-28,  1998. 


G.  WEB  SITE  INFORMATION 

A  web  site  is  currently  under  development.  It  will  contain  (i)  this  report,  (ii)  our  Matlab- based  virtual  analyzer 
application,  and  (iii)  links  to  sites  such  as  that  relate^  tp  Jhe  helicopter  data  [WOR], 


